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We present a recent analysis of e'/e in the 1/A'^c expansion. We show that the 1/Nc corrections to the matrix 
element of Qg are large and positive, indicating a A/ = 1/2 enhancement similar to the one of Qi and Q2 which 
dominate the CP conserving amplitude. This enhances the CP ratio and can bring the standard model prediction 
close to the measured value for central values of the parameters. Several comments on the theoretical status of 
e'/e and the errors in its calculation are given. 



1. INTRODUCTION 

Direct CP violation in K ^ tttt decays was 
recently observed by the KTeV and NA48 collab- 
orations pi, 0]. The present world average |0| for 
the parameter e'/e is Re e'/e = (19.3±2.4)-10-*. 
In the standard model CP violation originates in 
the CKM phase, and direct CP violation is gov- 
erned by loop diagrams of the penguin type. The 
main source of uncertainty in the calculation of 
e' /e is the QCD non-perturbative contribution re- 
lated to the hadronic nature of the K ^ tttt de- 
cay. Using the AS" — 1 effective hamiltonian. 
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u^C^{^i)Q^{^i) 



{iJL < rric), (1) 



the non-perturbative contribution, contained in 
the hadronic matrix elements of the four-quark 
operators Qi, can be separated from the pertur- 
bative Wilson coefficients Ci {fi) = Zi {p) + ryi (/i) 
(with r — —Xt/Xu and A, — V*g V_^). Introduc- 
ing {Qi)i = {{TT'K)i\Qi\K), the CP ratio can be 
written as 

s' Gf ujlraXt r I v^ ,^ , 

|_ 1 2^ y* {Qi)o 



2 |e|Re^o 



1-a 



(2) 



uj =ReAo/ReA2 = 22.2 is the ratio of the CP 
conserving K ^ ir-n isospin amplitudes; Q.i^ pa- 
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rameterizes isospin breaking corrections [g|. e'/e 
is dominated by {Qqjo and {Qs)2 which cannot be 
fixed from the CP conserving data Q, || . Beside 
the theoretical uncertainties coming from the cal- 
culation of the {Qi)i and of fJie, the analysis of 
the CP ratio suffers from the uncertainties on the 
values of various input parameters, in particular 
of the CKM phase in ImAt, of Age 
of the strange quark mass. 



A^, 



and 



2. ON THE COUNTING IN l/N^ AND 
THE USE OF THE LARGE-ATc VAL- 
UES FOR THE MATRIX ELEMENTS 

To calculate the hadronic matrix elements we 
start from the effective chiral lagrangian for pseu- 
doscalar mesons which involves an expansion in 
momenta where terms up to 0{p'^) are included 
Ig]. Keeping only terms of 0{p^) which con- 
tribute, at the order we calculate, to the K — > tttt 
amplitudes, for the lagrangian we obtain: 



C. 



eff = 

ll{{D^U^D>^U) + 4^(lnC/t - lnC/)2 

+ {XU^ + Ux^)) + Lsix^Ux^U + XU^XU^) 

+L,{D^U^D'^Uix^U + Uh)), (3) 

with (A) denoting the trace oi A, a — mi + rn^, — 
2w|., X = fM, and M = d\a.g{mu,md,ms). f 
and r are parameters related to the pion decay 
constant F-rr and to the quark condensate, with 



r = ~2{qq) / f'^ . The complex matrix [/ is a non- 
linear representation of the pseudoscalar meson 
nonet. The conventions and definitions we use are 
the same as those in 0, Jl Q]. The method we use 
is the \/Nc expansion |1C(| . In this approach, we 
expand the matrix elements in powers of the mo- 
menta and of l/Nc- For the l/N^ corrections we 
calculated chiral loops as described in refs. [|[ |9). 
Especially important to this analysis are the non- 
factorizable corrections, which are UV divergent 
and must be matched to the short-distance part. 
They are regularized by a finite cutoff Ac which 
is identified with the short-distance renormaliza- 
tion scale. The definition of the momenta in the 
loop diagrams, which are not momentum trans- 
lation invariant, is discussed in detail in ref. ||]. 
Other recent work on matrix elements in the l/N^ 
approach can be found in refs. [|ll], |l^, p^ . 

For the Wilson coefficients we use the leading 
logarithmic and the next-to-leading logarithmic 
values B. The absence of any reference to the 
renormalization scheme in the low-energy calcula- 
tion, at this stage, prevents a complete matching 
at the next-to- leading order |1J]. Nevertheless, 
a comparison of the numerical results obtained 
from the LO and NLO coefficients is useful as re- 
gards estimating the uncertainties and testing the 
validity of perturbation theory. 

As it is well known the large- TVc approximation 
fails completely in explaining the A/ = 1/2 rule. 
By taking the large-A^c values for the two domi- 
nant operators Qi and Q2 one obtains a A/ = 1/2 
CP-conserving amplitude which underestimates 
the data by roughly a factor of three. However, 
as pointed out in ref. [Eol, Q1.2 show an impor- 
tant specificity; they are expected to be largely af- 
fected by non-chirally suppressed corrections be- 
yond the large- ATc limit. In the counting in p^ 
and 1/Nc, this property is attributed to the fact 
that the tree level mesonic representation of Qi^2 
from the leading 0{p'^) chiral lagrangian intro- 
duces 0{p'^) terms with two derivatives. As a re- 
sult the one-loop contributions over these terms, 
which are 0{p^/Nc), are quadratically divergent, 
i.e. possibly very large because these corrections 
are not protected by any symmetry like the chiral 
one. Calculating the terms of 0{p'^/Nc) one ob- 
tains a large enhancement of the A/ = 1/2 ampli- 



tude |10(] in the range required to reproduce the 
experiment ||9(]. This illustrates how important 
are the non-chirally suppressed 1/A'c corrections 
for an understanding of the K — )■ tttt amplitudes. 
In this context, for e'/e, it is very important to 
investigate the counting in p^ and 1/iVc for the 
dominant operators Qe,8 and to compare it with 
the one for Qi^2- For Qg the mesonic representa- 
tion from the 0{p^) lagrangian is 0{p'^) and does 
not have any derivative (because Qs is a density- 
density operator whereas Qi,2 are of the current- 
current type). As a result we do not expect 
any quadratic term but only chirally suppressed 
logarithmic or finite terms. The large-A'c limit, 
Sg = 1, is therefore expected to be a much 
better approximation than for Qi,2- Different is 
the case of Qq due to the fact that there is no 
tree level contribution from the leading 0{p^) la- 



grangian, and the large- iVc value B^ 
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to the tree level contribution from the 0{p'^) la- 
grangian which, as for Qi^2, carries two deriva- 
tives [i.e. is C(p^)]. Therefore, as for Qi,2, we ex- 
pect the large- iVc value of Sg ' ' to be affected by 
large O^p^/N^) corrections [- 0(100%)] result- 
ing from the loops over the 0{p'^) tree operator. 
This shows clearly that, in the same way as for the 
A/ = 1/2 rule, no clear statement can be done 
on the expected size of e'/e without calculating 
these 1/A^c non-factorizable corrections. 



3. ANALYSIS OF e'/e 

Analytical formulas for all matrix elements, at 
next-to-leading order in the twofold expansion in 
powers of momenta and of 1/Nc, are given in 
refs. [pi P[. In the pseudoscalar approximation, 
the matching has to be done below 1 GeV. Vary- 
ing Ac between 600 and 900 MeV, the bag factors 
BJ^/^^ and S^^^^^ take the values 8.2 - 14.2 and 
2.9 — 4.6; quadratic terms in (Qi)o and {Q2)o 
produce a large enhancement which brings the 
A/ =1/2 amplitude in agreement with the data 
[^. Corrections beyond the chiral limit were 
found to be small. 

For {Qe)o and {Qs)2 the leading non- 
factorizable loop corrections, which are of 
0{p^/Nc), are only logarithmically divergent |a|. 



Including terms of 0(p"), 0(p2), and 0{p°/Nc), 
B^^/^^ and B^^/^^ take the values 1.10 - 0.72 and 
0.64 — 0.42. As a result, as with large-TVc val- 
ues, e'/e is obtained generally much smaller than 
the data, except for input parameters taken close 
to the extreme of the ranges we considered. How- 
ever, as stated above, since the leading 0{p^) con- 
tribution vanishes for Qq, corrections from higher 
order terms beyond the 0{p^) and 0{p^ /Nc) are 
expected to be large. In ref. 0| we investigated 
the 0{p'^/Nc) contribution, i.e., the 1/iVc cor- 
rection at the next order in the chiral expan- 
sion, because it brings about, for the first time, 
quadratic corrections on the cutoff. From count- 
ing arguments and more generally from the fact 
that the chiral limit is assumed to be reliable, 
the quadratic terms (which are not chirally sup- 
pressed) are expected to be dominant. It is still 
desirable to check that explicitly by calculating 
the corrections beyond the chiral limit, from log- 
arithms and finite terms, as done for Qi and Q2. 
Numerically, we observe a large positive correc- 
tion from the quadratic term in (Qelo- This point 
was already emphasized in ref. |1^. The slope 
of the correction is qualitatively consistent and 
welcome since it compensates for the logarithmic 
decrease at 0{jP /Nc). Varying Ac between 600 
and 900 MeV, the B^ factor takes the values 
1.50 — 1.62. Qq is a AJ — 1/2 operator, and 
the enhancement of {Qe)o indicates that at the 
level of the 1/Nc corrections the dynamics of the 
A/ = 1/2 rule applies to Qq as to Qi and Q2. 
The size of the enhancement for Bq appears 
however to be smaller for Qq than for Qi^2 due to 
a smaller coefficient of the quadratic term. This 
coefRcient is nevertheless large enough to increase 
e' /e by almost a factor of two. 

Using the quoted values for Bq ' ' together 
with the full leading plus next-to-leading order B 
factors for the remaining operators J 3] the results 
we obtain for e'/e are given in tab. |l for the three 
sets of Wilson coefficients LO, NDR, and HV and 
for Ac between 600 and 900 MeV. The numbers 
are close to the measured value for central values 
of the parameters (upper set) . They are obtained 
by assuming zero phases from final state inter- 
actions. This approximation is very close to the 



Table 1 

Numerical values for e'/e (in units of 10~*) as 

explained in the text. 



LO 


14.8 < e'/e < 19.4 




NDR 


12.5 < e'/e < 18.3 


central 


HV 


7.0 < e'/e < 14.9 




LO 


6.1 < e'/e < 48.5 




NDR 


5.2 < e'/e < 49.8 


scanning 


HV 


2.2 < e'/e < 38.5 





results we would get if we used the small imagi- 
nary part obtained at the one-loop level 0. 

Performing a scanning of the parameters 
[125 MeV < 771^(1 GeV) < 175 MeV, 0.15 < 
^B < 0.35, 1.04 • 10-4 < ImAt < 1.63 • 10-^ 
and 245 MeV < Aqcd < 405 MeV] we obtain the 
numbers in lower set of tab. y. They can be com- 
pared with the results of refs. |, ^ |l^, [l|, |l8) . 
The values of Bg can also be compared with 
refs. |1^, [ij]. Other recent calculations are re- 
ported in refs. |g^, |l| [2^. The wide ranges re- 
ported in the table can be traced back, to a large 
extent, to the large ranges of the input parame- 
ters. This can be seen by comparing them with 
the relatively narrow ranges obtained for central 
values of the parameters. The parameters, to a 
large extent, act multiplicatively, and the large 
range for e'/e is due to the fact that the cen- 
tral value(s) for the ratio are enhanced roughly 
by a factor of two compared to the results ob- 
tained with B factors for Qg and Qg close to the 
VSA. More accurate information on the param- 
eters, from theory and experiment, will restrict 
the values for e'/e. 

To estimate the uncertainties due to higher 
order final state interactions we also calculated 
e'/e using the real part of the matrix elements 
and the phenomenological values of the phases 
11,(50 = (34.2 ± 2.2)° and ^2 = (-6.9 ± 
0.2)°, i.e., we replaced |E^y.(Q^>/| in Eq. §) 
by ^^yiRe((5i)7/coS(y/. The corresponding re- 
sults are given in tab. 0. They are enhanced by 



Table 2 

Same as in Tab. H, but now with the phenomeno- 

logical values for the phases. 



LO 


19.5 < e'le < 24.7 




NDR 


16.1 < e'/e < 23.4 


central 


HV 


9.3 < e'/e < 19.3 




LO 


8.0 < e'/e < 62.1 




NDR 


6.8 < e'/e < 63.9 


scanning 


HV 


2.8 < e'/e < 49.8 





~ 25 % compared to the numbers in tab. [I[ To re- 
duce the FSI uncertainties in the l/Nc approach 
it would be interesting to investigate the two-loop 
imaginary part. By doing so we expect to get 
phases very close to the ones of ref. |24 which 
have been obtained in Chiral Perturbation The- 
ory at the same order and reproduce relatively 
well the data. In this sense we expect to get re- 
sults close to the ones of tab. 0. A comparison 
of tabs. and will be however still useful to es- 
timate higher order corrections (e.g. for the real 
part). In our analysis part of the uncertainty from 
higher order corrections is also included in the 
range due to the (moderate) residual dependence 
on the matching scale. In order to reduce the 
scheme dependence in the result, appropriate sub- 
tractions would be necessary [^ |l^, |5| . Finally, 
it is reasonable to assume that the effect of the 
pseudoscalar mesons is the most important one. 
Nevertheless, the incorporation of vector mesons 
and higher resonances would be desirable in or- 
der to improve the treatment of the intermediate 
region around the rho mass and to show explic- 
itly that the large enhancement we find at low 
energy at the level of the pseudoscalars remains 
up to the scale ^ rric, where the matching with 
the short-distance part can be done more safely. 

4. ON THE SIZE OF THE ERRORS IN 
THE ANALYSIS of e'/e 

As shown in tabs. ^ and the errors we obtain 
for e'/e are large. We believe however that the 



uncertainties are not largely overestimated and 
reflect well our present knowledge on e'/e. 

Note that we performed a scanning of the pa- 
rameters. Refs. g |l^, |l^, Q used in addition 
to the scanning method also a Monte Carlo pro- 
cedure with gaussian distributions for the exper- 
imental input and flat distributions for the theo- 
retical parameters. In this way a probability dis- 
tribution is obtained for e'/e, and the authors 
gave the median and the 68 % confidence level in- 
terval. We would like to emphasize that the use 
of this C.L. interval removes part of the hadronic 
uncertainties and leads to a range for e'/e which 
is two times smaller than the one obtained from 
a full scanning over the ranges of the theoreti- 
cal (mostly hadronic) parameters. In our opinion 
the Monte Carlo analysis is misleading because 
there is no justification for assuming any proba- 
bility distribution for theoretical parameters, and 
it leads to an underestimate of the uncertainties 
in the calculation. To illustrate this point one 
might note that values for e'/e above the 68% 
C.L. range given e.g. in ref. iQ can be obtained 
for central values of the experimental parameters 
and for quite reasonable values of the theoretical 
ones within the ranges considered in this refer- 
ence. Therefore we think that a full scanning of 
the theoretical parameters, with gaussian distri- 
bution or full scanning for the experimental pa- 
rametersg, gives a better idea of the uncertainties 
in the CP ratio. A similar comment applies to 
the range for ImAt since the determination of the 
CKM phase involves many non-perturbative the- 
oretical parameters (as e.g. Bk)- 

Moreover, as we explained above, non-chirally 
suppressed corrections beyond the large Nc limit 
Therefore, in our opinion 



(1/2) 
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are essential for B, 

larger errors for Bq'''^' should have been taken in 
ref. |1^ since the authors did not calculate this 
parameter (in this case, from the counting in p^ 
and l/Nc, the errors could be taken as large as 
100%). The same comment applies to ref. p^ 
where a similar range for {Qe)o was adopted 
[i?g ' ' was varied around its large- A^c value tak- 



^ Since the ranges in tabs, hi and Marc predominantly due to 
the theoretical uncertainties, the difference between these 
two procedures is moderate. 



ing an error of 100 %, but rris was fixed adopting 
the value (m^ + md)(2 GeV) = 130 MeV without 
considering any error on it]. A similar statement 
applies to the dispersive analysis of ref. [gg which 
does not give access to any of the scale depen- 
dent non- fact or iz able terms (the non-chirally sup- 
pressed quadratic corrections in particular). The 
calculation of the scale dependent terms is not 
easier for lower values of the squared momentum 
of the kaon, and dispersive techniques do not help 
in their calculation. To neglect these terms leads 
to a failure in reproducing the A/ = 1/2 rule and 
is also not justified for Qe-QJ As pointed out in 
ref. fill, in the Chiral Quark Model the corre- 
lation between the A/ =1/2 amplitude and e'/e 
(used to fix the parameters necessary to estimate 
the matrix element of Qq) is subject to poten- 
tially large uncertainties. As for flj^-, a minimum 
error of '^ 0.10 seems to be required for a careful 
estimate of e'/e. The value fiig = 0.16 ± 0.03 in 
ref. 1^ (which has been used in ref. Q) was ob- 
tained by investigating the it — ij contribution to 
this parameter (including rj — -q' mixing). How- 
ever, corrections beyond this term could be non- 
negligible as suggested by the numerical results 
of refs. H, H. 

We would like to emphasize also that the ~ 
25 % error obtained by comparing tabs. |] and Q 
should be included by any analysis which either 
does not include final state interactions or does 
not reproduce well the numerical values of the 
phases. This estimate of part of the neglected 
higher order corrections is usually not taken into 
account. 

We conclude that at present there is no method 
which can predict e'/e with an error much smaller 
than the one presented in tabs. and which 
give a good idea of the uncertainties involved in 
the calculation of the CP ratio. The statement, 
that the experimental data can be accommodated 
(only) if all the hadronic parameters are taking 
values at the extreme of their reasonable ranges 
(see e.g. ref. [pOl), which is based on the use of 
the 68% C.L. intervals of refs. |, |6[ |l|, has 
only weak theoretical foundations. 



^Further comments on the dispersive analysis of the FS I 
effects in the calculation of e'/e can be found in ref. [G7| . 



5. COMMENTS ON e'/e 

Following our analysis, described above, a se- 
ries of comments can be made: 

• The use of the large- iVc value or values close 
to it is not justified for Qq in e'/e in the 
same way as for Qi,2 in the A/ = 1/2 rule. 

• The related claim that we expect in gen- 
eral in the standard model a value of e'/e 
smaller than the data is therefore not justi- 
fied. One should note that the main meth- 
ods used to calculate the 1/Nc corrections 
H, 111, 11^ all find them large and positive. 

• Dispersive techniques as proposed in 
ref. pq] do not help in the calculation of 
many of these corrections (i.e. of the scale 
dependent terms). 

• Errors are large. A value between ~ few • 
10^"* and ~ 5 • 10"'^ appears perfectly plau- 
sible without going beyond the reasonable 
ranges of the parameters. 

• Therefore e'/e cannot be used at present 
to investigate new physics. Even if ImAt 
was found as small as 0.6 • 10^^ as could 
be suggested by recent measurements [ pi[ , 
we could still not exclude that the standard 
model reproduces the data. 

• At low energy there is a significant enhance- 
ment of -Bg [R|. It is quite reasonable 
that below 500 - 600 MeV the model inde- 
pendent lagrangian of Eq. (^ gives the bulk 
of the result. This is an important indica- 
tion for a large value of e'/e in accordance 
with the data. The results of refs. ||l|, |l7) 
also point towards this direction. 

• The effects of dimension-eight operators 
could possibly change the results largely 
P^ and should certainly be investigated. 

• Despite the recent progress in the calcula- 
tion of 17iB [pi pSl p^ the problem of an ac- 
curate calculation of this parameter is still 



relevant in the same way as for B^ 



(1/2) 



Questions (E. Pallante, Univ. Barcelona): 



1) Can you comment on the one-loop estimate of 



Bf'^'^ in ref. ll^? Tlie results of ref. fig (which 
have been obtained in the chiral limit) and our re- 

(3/21 

suits for B^ are not incompatible within their 
respective errors especially for moderate values 
of Ac. 

2) The value of r2iB = 0.16 ± 0.03 is actually 
not a pure large-Nc calculation. There is no rea- 
son to expect large 1/Nc corrections in this case. 
There is a priori no reason for the 1/Nc correc- 
tions from irreducible one-loop diagrams beyond 
the reducible ones calculated in ref. H to be neg- 
ligible (see also refs. |8[ |9|). 
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